Multipartite entanglement generation and fidelity decay in disordered qubit systems 
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We investigate multipartite entanglement dynamics in disordered spin-1/2 lattice models exhibit- 
ing a transition from integrability to quantum chaos. Borrowing from the generalized entanglement 
framework, we construct measures for correlations relative to arbitrary local and bi-local spin ob- 
servables, and show how they naturally signal the crossover between distinct dynamical regimes. 
Analytical estimates are obtained in the short- and long-time limits. Our results are in qualita- 
tive agreement with predictions from random matrix theory and relevant to both condensed-matter 
physics and to the stability of quantum information in disordered quantum computing hardware. 
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Developing a quantitative understanding of the struc- 
tural and dynamical properties of entanglement in many- 
body quantum systems is a critical challenge for both 
condensed-matter theory and quantum information sci- 
ence. Low-dimensional disordered spin models offer, in 
this context, an ideal testbed for theoretical analysis. On 
one hand, these systems are simple enough for analytic 
benchmarks to exist in limiting situations, yet capable to 
demonstrate a broad typology of complex quantum phe- 
nomena. The latter range from field- or disorder-driven 
structural ground-state changes responsible for quantum 
phase transitions Q, to dynamical crossovers from in- 
tegrable to non-integrable regimes and the emergence of 
quantum chaos 2\. On the other hand, arrays of interact- 
ing spin- 1 /2s naturally schematize quantum computing 
hardware as diverse as semiconductor quantum dots [3|, 
superconducting Josephson qubits Q, electron floating 
on He |5| , and optical lattices |6| - disorder resulting from 
the unavoidable presence of imperfections in the single- 
qubit energy spacings and/or inter-qubit couplings. 

Following 0, intense effort has been devoted to both 
assess the impact of disorder on quantum computing per- 
formance 0, 0, [ljj, an d to characterize entanglement 
across the ensuing transition to quantum chaos [111 |l2| • 
These studies point to the key role of spectral properties, 
as captured by the Local Density of States (LDOS) [l3| . 
in determining the system stability against the disorder. 
Changes in the LDOS profile are ultimately responsible 
for the existence of different dynamical regimes, as re- 
flected by qualitative changes in the rate of fidelity de- 
cay [l II . Numerical experiments confirm that the same 
regimes translate into distinctive static 01 an d dynamic 
|l5| properties of pairwise entanglement, as quantified by 
concurrence |l6). While providing suggestive evidence, 
such analyses are not fully satisfactory for two reasons. 
First, concurrence lacks a direct physical interpretation, 
hindering the possibility to relate entanglement to spec- 
tral properties. Second, concurrence is a bipartite mea- 
sure, preventing the quantification of multipartite corre- 
lations which dominate in strongly coupled scenarios. 



In this Letter, we overcome the above limitations by 
exploiting Generalized Entanglement (GE) [l7| as a set- 
ting for defining entanglement relative to arbitrary sets 
of observables. In (is| . GE measures constructed from 
algebras of fermionic operators have been applied to the 
study of broken-symmetry quantum phase transitions in 
exactly solvable models, notably the spin-1/2 XY chain 
in a transverse field. Here, we show how GE contributes 
to the understanding of standard multipartite correla- 
tions between distinguishable systems, by selecting local 
and bi-local algebras of observables corresponding to sin- 
gle and pairs of qubits, respectively. Beside providing a 
transparent interpretation of the results based on con- 
currence [l2l llfij . our approach allows for quantitative 
insight about entanglement dynamics starting from arbi- 
trarily correlated initial states. We quantify the influence 
of LDOS properties in two limits: at short-time, by es- 
tablishing a direct link with fidelity decay; at long times, 
by relating the entanglement saturation value to the par- 
ticipation ratio of the asymptotic many-body state. 

The model- We focus on a two-dimensional lattice of n 
disordered spin-1/2 particles (n even) described, in units 
H = 1, by the following Hamiltonian: 



H 



E 

i=i 



A + 6j]a 



2> 



H A + H S + Hj 



where Oa\ a £ {0, x, y, z}, <7q 1> — I, are Pauli operators, 
and (i,j) stands for nearest-neighbor sites. Open bound- 
ary conditions and energy units A = 1 are used. The 
parameters Sj , Jy characterize the disorder in the on-site 
energy splittings and two-body couplings, respectively. 
We assume that Sj, Jij are uniformly random over [—6, S], 
[—J, J], with 5, J > 0. The above model, which belongs 
to the random transverse-field Ising lattice class 0] , was 
recently used to describe a quantum register with static 
imperfections • Understanding the influence of the lat- 
ter may be especially important for controlling entan- 
glement in quantum computers with always-on interac- 
tions |2fl | and one-way quantum architectures [2ll |. 
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FIG. 1: (Color Online) Local purity critical time t c (Pi(t c ) — 
0.9) as a function of J/5 for 6 = 0.02, 0.05, 0.1, 0.2 (diamonds, 
triangles, circles, squares). To smooth statistical fluctuations, 
here and in the following we consider averages {Pi(t))D over a 
number N r = 10 of disorder realizations. We assume n = 10 
unless otherwise specified. Inset: Time evolution of Pi(t) 
in the FGR J = 8/10 = 0.01 (red) and ergodic J = 6 = 
0.1 (blue) regime starting from the state of Eq. © for n = 
14, n_B = 0. Dotted lines: Exponential and Gaussian fits, as 
in Eq. ©. 



The spectrum of is composed of n + 1 levels, 
with energy Eb — A(2fc — n) and degeneracy Ns(k) = 
n\/[k\(n - k)\], k = 0,...,n. When 5, J -C A, the degen- 
eration is removed and n + 1 bands appears. The width 
A B of the bands depends on the ratio J/5 Q • In the limit 
where J <C 6, A s ~ as the spread is led by the diag- 
onal term H$. Correspondingly, bands are Gaussian. In- 
side each band, the effective number of states coupled by 
Hj is given by the width of the LDOS. For J < 5/n = J c , 
Hj couples few states, and the system is slightly per- 
turbed, while for J > J c quantum chaos sets in and 
the unperturbed levels are coupled to a quasi-continuum 
set of states. The LDOS becomes Lorentzian, with a 
width determined by the Fermi Golden Rule (FGR), 
r F = J 2 n/8 0. Increasing J, when Tp ~ Ag, that 
is for J ~ Je — S/n 1 / 4 , all the levels inside a given 
band are mixed and the LDOS approaches the level den- 
sity. Thus, the LDOS becomes a Gaussian with width 
~ J pj. In summary, three distinct regimes exist 
- perturbative (J < J c ), FGR (J < Je), and ergodic 
(Je < J < A). It has been shown that starting from an 
eigenstate \tpo) of Ha, the survival probability (fidelity 
henceforth, as for initial eigenstates the two quantities 
coincide) F(t) — | (^olV't) | 2 also follows three different 
behaviors: F(t) oscillates near one in the perturbative 
regime, while it decays as an exponential or a Gaussian 
in the FGR and ergodic regime, respectively Q. The 
connection between LDOS shape and fidelity decay has 
been unveiled in |8( by means of a perturbative relation 
between F(t) and the Fourier transform of the LDOS. 



Local and bi-local purities.- An initial pure state |^o) 
evolving under H remains pure for any fixed disorder 
realization. Accordingly, we consider pure-state entan- 
glement throughout. The basic intuition underlying the 
GE notion pjj is to quantify how entangled a state \ip) is 
relative to an observable set O in terms of how pure \ip) 
remains upon restricting operational access to O. While 
no single set O can exhaust the complexity of multipar- 
tite entanglement, two simple choices will illustrate the 
usefulness and flexibility of the algebraic approach. First, 
we consider the set O = 0\ generated by arbitrary local 
observables that is, we probe the average entanglement of 
each qubit with the rest of the lattice. The corresponding 
(time-evolved) local purity measure is defined as 

i— l,ri 

W» = ~ E K^hMI 2 . (i) 

a— x,y,2 

For this choice of observables, GE coincides with global 
multipartite entanglement as quantified by the Meyer- 
Wallach metric Q |30> with <3 = 1- Pi0& 
As a second observable set, whose physical motivation 
will become clear later, we choose O = O2 generated by 
all the observables acting on pairs of nearest neighbor 
spins. We compute a bi-local purity measure as 

where the tilde means that the identity term with a = 
(3 = is omitted. P2 may be thought as a "coarse- 
grained" version of P±, resulting from ignoring the fine 
structure given by arbitrary correlations within each pair. 
Both measures are normalized to give one (zero) on states 
which arc fully separable (contain maximal GE) relative 
to the corresponding algebra. Remarkably, both Pi and 
Pi are directly measurable quantities in principle |23| . 

Initial separable state - We first study multipartite en- 
tanglement generation starting from a separable state 
in the central band (k = n/2) that is, |V>o) = 
1010101 . . .01) = |c), where {0, 1} label the states of each 
spin in the computational basis and c is the integer given 
by the corresponding binary string. The evolution of 
Pi (t) under different values of J at fixed 6 are depicted in 
the inset of Fig. ^ The decay of Pi (t) clearly follows two 
different dynamical laws. Specifically, data in the FGR 
and in the ergodic regime are described, respectively, by 

P[ (t) » e- CI > 4 , if (t) « e - c " r - t2 , (3) 

where the constants C, C depend on the initial state and 
the lattice coordination number (see e.g. [j^|). Let t c be 
the time it takes for Pi (or P2) to reach the value K e.g., 
Pi(t c ) = K: We find that t c F - 1/J 2 in the FGR regime, 
whereas t E ~ 1/J. These behaviors have been verified 
over a wide range of disorder parameters, see Fig. ^ 
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FIG. 2: (Color Online) Saturation values of 1-1 (empty sym- 
bols) and of (Pf 1 -1)-R (filled symbols) vs J/8 for 5 = 0.05 
(circles), 8 = 0.1 (squares), 5 = 0.2 (diamonds), and rescal- 
ing factor R — 7. The straight line is proportional to 
C{J/5) 2 N B (k = n/2), C ~ 0.7. Inset: Saturation value of 
the IPR (triangles up) and Pj -1 (triangles down) in the er- 
godic regime (8 = 0, J = 0.1) as a function of n. The dashed 
line gives Nb(ti/2) vs n. 



A simple physical interpretation of the above results 
follows from the possibility to directly relate GE dynam- 
ics to fidelity decay for any computational eigenstate |c). 
While a full derivation will be presented elsewhere, the 
key steps are (i) to realize that only o* observables con- 
tribute to the evolution of Pi (t) (ii) to isolate the fidelity 



term |(c|e" 



-iHt\ 



in the resulting z-purity. This yields 



1 n 



3 = 1 



where ay(i) = £ p/c |<p|e- ,Ht |c)| a (-) ,r 'to and {~Y liq) 
= ( — ) L<?/2 J J i s the parity of the jth qubit in the com- 
putational state \q). One can show that each term ctj(t) 
is of order 0((J/S) 2 t 2 ). Thus, according to the above 
equation, the connection between local purity and fidelity 
decay (hence LDOS via the relation found in ||) becomes 
exact in the limit t — > 0. For sufficiently short times, the 
first term still dominates and the dynamics is governed 
by F(t) 2 , whereby the two regimes of Eq. © arise. 

For times much longer than r^^, the asymp- 
totic GE amount may be estimated using random 
matrix theory J2j. By working within the Gaus- 
sian orthogonal ensemble, we represent the many-body 
state IV'oo) as a random superposition of Noo unper- 
turbed states within the central band that is, |V>oo) = 
^2p=iW p \p), with amplitudes w p satisfying the nor- 
malization constraint. For sufficiently large A^, cor- 
relations between different components may be ne- 
glected, and the probability density of each component 
is well approximated by the Porter-Thomas distribu- 



tion, P PT (K| 2 ) = {2TTN 00 \w p \ 2 )- 1 / 2 eM-N 



72). 




FIG. 3: (Color Online) Local purity Pi(i) versus time in the 
FGR (black lines) and ergodic regime (red lines) for different 
initial states with ub = 0,1,2,3, 8 = 0.1, and J = 0.01 
(black curves), J = 0.1 (red curves). Dashed lines: Different 
separable initial states \ip ) = |0101010110), |010101010) ® 
(|0) + |1})/V2 for 8 = 0.1, J = 0.01. Inset: The sames curves 
rescaled by Pi(0). Also shown (dot-dashed lines) is the local 
purity decay in the ergodic and FGR regimes for an initial |W) 
state. 



yields (P{°)d = <-Pi°°)pt = consistent with the 

scaling 3/(A^oo + 1) predicted for a random pure state 
under the invariant Haar measure |25|. 

Let the degree of derealization of IV'oo) be quantified 
by the Inverse Participation Ratio (IPR) 0, 0, E] that 
is, 1= l/£, where the asymptotic participartion num- 
ber = £ P KI 4 - By estimating = 3/AU 

(P^°) D = 2/3(£°°)d - thus, GE directly quantifies de- 
localization in this regime. Similarly, (I°°)d ~ 

— 1/2 

up to deviations of the order 0(Noq ), consistent with 



Gaussian fluctuation statistics |2|. Because the dynam- 
ically accessible states are determined by the ratio be- 
tween the LDOS width and the bandwidth, we further 
estimate as (Tp/ Ab)Nb for any initial eigenstate in 
the central band and in the FGR regime, whereas the 
IPR is constant in the ergodic regime. These predictions 
have been confirmed numerically (see Fig. [2J ^26] . Thus, 
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A, 



'J )d T F N B 



(4) 



Evaluating the disorder-averaged asymptotic local purity 



As shown in Fig. El Eq. © qualitatively agrees with data 
over an extensive range of parameters. Note that these 
asymptotic results, as well as the initial decay laws of Eq. 
(|3*|l. are valid for the dynamics of any initial separable 
state, not just of a computational state (see Fig.|3). 

Initial entangled states - The analysis may be extended 
to the dynamics of arbitrarily entangled initial states. 
Consider a state in the central band containing only bi- 
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partite entanglement first, for instance 



n — 2rif 



\Mn B )) = |01...01>® 



T2 m 



|io» 



(5) 



where rig is the number of Bell pairs, ng = recovering 
the fully separable case. The local purity evolution is de- 
picted in Fig.|3| Reflecting the fact that Pi is sensitive to 
all correlations between spins, the initial value is lower 
the larger ns, Pi(0) = 1 — 2ns /n, using Eq. 0. For 
t > 0, Px(t) decays similarly to the separable case, two 
distinct dynamical regimes emerging for jib — 0,1,2,3. 
As expected, data corresponding to a given regime ap- 
proximately fall on the same curve once rescaled by Pi (0) 
(Fig- El inset). This may be quantitatively understood 
by studying the evolution of the bi-local purity defined 
in Eq. PJl. By construction, Pi is insensitive to any 
pairwise correlation present in the state (J5J, effectively 
mapping the analysis back to the separable case ns = 0. 
P 2 (t) is depicted in the inset of Fig. 01 Clearly, P 2 (0) = 1 
irrespective of ub- The time decay then follows the two 
regimes predicted by ©, again reflecting the underlying 
structure of the LDOS (see Fig. 0J. 

As a representative example of an initial state con- 
taining genuine multipartite entanglement, we focus on 
a so-called |W) state that is, an equally weighted super- 
position with equal phases of the Nb(1) — n states with 
^-magnetization J2j a< z = n — 2. Within our disorder 
model, |W) is automatically protected against the effects 
of H$, Hj as long as the coupling between different bands 
remains small, as assumed so far. Indeed, the local pu- 
rity remains almost constant (data not shown). Thus, 
to analyze how the (tripartite) correlations contained in 
a |W) are affected by static random imperfections we set 
A = 27]. Again, the same qualitative picture arises, 
with two distinct dynamical regimes (inset of Fig. an d 
two different scalings of critical decay times t c (Fig.0J. 

Conclusions - We established the existence of distinc- 
tive signatures in the dynamics of multipartite entangle- 
ment in a quantum many-body system subject to static 
disorder. While a more systematic analysis is certainly 
desirable, we expect our main conclusions to prove valid 
under more general conditions, including different dimen- 
sionality and/or model Hamiltonians, as well as higher- 
spin systems. Beside reinforcing the usefulness of the GE 
notion as a diagnostic framework for complex quantum 
systems, the deep connections between GE, LDOS, and 
fidelity decay emerging from our work are likely to have 
broader implications across the fields of quantum infor- 
mation, condensed-matter physics, and quantum chaos, 
stimulating crosstalk between different communities. 

We thank W. G. Brown, R. Fazio, G. Ortiz, L. F. 
Santos, and Y. S. Weinstein for feedback. Partial sup- 
port from the IST-SQUBIT2, from IBM (Faculty Awards 
2005), and from Constance and Walter Burke through 
their Special Projects Fund in QIS is acknowledged. 
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FIG. 4: (Color Online) Purity critical time t c as a function 
of J/8 for different initial states, ub = 1 (circles), ub = 2 
(squares), ub — 3 (diamonds), |W) (blue triangles,) and dif- 
ferent purity measures: Pi(t) (empty) and Pz(t) (full). In- 
set: Piit) vs time in the FGR regime (black lines) and in 
the ergodic regime (red lines) for different initial states with 
n B =0,1,2,3. 
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